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Abstract 

A novel method to calculate fs on the lattice is introduced, based on the study of the dependence of 
finite size effects upon the heavy quark mass of flavoured mesons and on a non-perturbative recursive 
finite size technique. We avoid the systematic errors related to extrapolations from the static limit or 
to the tuning of the coefficients of effective Lagrangian and the results admit an extrapolation to the 
continuum limit. We perform a first estimate at finite lattice spacing, but close to the continuum limit, 
giving f B = 170(11)(5)(22) MeV. We also obtain f Bg = 192(9)(5)(24)MeV. The first error is statistical, 
the second is our estimate of the systematic error from the method and the third the systematic error 
from the specific approximations adopted in this first exploratory calculation. The method can be 
generalized to two-scale problems in lattice QCD. 
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1. Introduction 



Lattice QCD evaluations of quantities characterised by two scales with a large hierarchy require 
in general a very high lattice resolution and a sizeable total physical volume to correctly account 
the dynamics of the small distance scale and to dispose of the finite size effects related to the large 
distance scale. A good example is provided by the pseudoscalar B meson decay constant [1], 
where the small distance scale is represented by the inverse of the bottom quark mass and the 
large distance scale by the radius of the B meson, related in turn to the inverse of the light quark 
mass. A straight evaluation of the decay constant would require lattices with N = 80 4 points or 
more, exceeding the present generation computers capabilities, and, in the case of unquenched 
simulations, the ones of the next generation. One resorts to approximate calculations based on 
extrapolations from the static limit or on non-relativistic formulations of standard QCD. All 
the available methods introduce systematic errors related to extrapolation fits and/or to the use 
of effective Lagrangians. We present a novel approach based on the study of the dependence 
upon the heavy quark mass of finite size effects for the pseudoscalar decay constant of heavy 
flavoured mesons. The basic assumption is that the finite size effects are mainly related to the 
light quark mass and rather insensitive to the one of a sufficiently heavy quark. We discuss 
the general features of the method assuming the continuum limit has been taken. Corrections 
specific to the finite lattice spacing calculation presented in this first paper are discussed later. 
The relevant quantity is the ratio a of the pseudoscalar constants at different volumes: 

_ f B (2L) 

a = 1WY (L1) 

where /b(L) is the value of the decay constant on a volume with linear size L. The dimensionless 
a depends on general grounds upon three dimensionless variables: m^L, m h L and Aq CD L. For 
a sufficiently large heavy quark mass m/,, the dependence is basically dominated by the light 
quark and the expansion for large takes the form 

a = a{miL, K QC dL) + — (1.2) 

A simple phenomenological ansatz for a can be made based on the concept of a reduced mass 
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constructed out of the heavy and light quark masses 



a = a 



(m red L, KqcdL) 



(1.3) 



where 



(1.4) 



faired 



Ml + M2 



The quantity /x^ is a function of the quark mass, but not only: indeed, for very light masses, 
finite size effects are regulated by the physical meson size, which is expected to remain finite 
when the light quark mass tends to zero. We will show later some evidence for ^ being a simple 
linear combination of the light quark mass and Aqcd- 

A crucial question is the threshold value of the quark mass on a given volume where the large 
rrih expansion becomes reliable. As we will see, this value falls in a mass range of the order of a 
couple of GeV in the renormalization invariant mass scheme, where the calculation on a single 
lattice is affordable. Under these circumstances, the strategy to obtain is the following. One 
first performs a calculation on a lattice where the resolution is suitable for b quark propagation, 
but the total volume is unavoidably a small one. This sets on a finite volume. In order to 
connect to the large volume results, one needs the step scaling function a for values of heavy 
quark masses generally lower than those of the simulation where the finite size value of fs was 
obtained. The possibility of extrapolating a to heavier masses depends upon the validity of the 
asymptotic expansion: in a favourable will be the real one, one can evaluate the finite 

size effects in a reliable way, connecting, by a repeated iteration of the procedure, small volume 
values of /b to the ones on large volumes, 



and the recursion stops on a volume where a ~ 1 within a required precision. The continuum 
limit is obtained by extrapolating to zero lattice spacing the step scaling function obtained at 
fixed physical quantities. This paper deals with a first exploration of the method at finite lattice 
spacing, suitably chosen to limit the systematic errors from lattice artifacts. Sec. 2 is dedicated 
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f p B hys = f B (L ) a(L ) a(2L ) ... 



(1.5) 



to the description of the general aspects of the calculation, sec. 3 to its specific details and to 
the results plus some comments. 



2. Theoretical framework 

This section is meant to set the notation, define the lattice observables and describe the strat- 
egy of the calculation. A recipe for defining the heavy-light states on a finite volume is then 
discussed. 

2.1 General strategy 

The calculation of fs is done on a set of lattices with topology T x L 3 in the Schrodinger 
Functional scheme [2,3], where gauge and fermion fields fulfill periodic boundary conditions 
along the space directions and Dirichlet boundary conditions at the beginning and at the end 
of the lattice history, and the following set of parameters is used 



Here C and C are the boundary gauge fields and 9 is a topological angle which enters into the 
definition of the Schrodinger Functional. In order to have a safe extrapolation to the continuum, 
non pertubatively 0(a) improved action [4] and operators are used. Within this framework, the 
gauge invariant correlation functions which describe the propagation of a heavy-light pseu- 
doscalar meson from the low boundary to the bulk and across the two boundaries are 



T = 



2L, 



C = C = 



9 = 



(2.1) 



U(x ) = -y<^(*)<S), h = -\(SS') 



(2.2) 



where the operators A { 



o(x), S and S' 



interpolate the meson field and are given by 




(x) = tp h (x)j^ 5 'tpe(x) + 



ac A 



(2.3) 



2 




(2.4) 
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We take the 0(a) improvement coefficient ca, which appears in the axial current, from [4]. The 
quantum mechanical representation of the correlations functions (2.2) can be found in [5]. Here 
we report for convenience their large time behavior, referring the reader to that paper for the 
notation: 



L 
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fA(x ) = —P (0,0\A \0,hi) exp(-x M M ) 



x 



1 + v B e - XO A +r] Ae 



-(T-xo)mc 



(2.5) 



fi = X -p 2 exp(-TM M ) (2.6) 

The state \0,h£) represents by definition the lowest heavy-light eigenstate with the quantum 
numbers of a pseudoscalar, and the matrix element (0, 0| Ao |0, hi) is related to the heavy-light 
meson decay constant through the relation 

Z A (0, 0|A)|0, hi) = f M M M {2M M L 3 )- 1 ' 2 (2.7) 

with the improved axial current renormalisation constant Za given by 

Z A = Z A (1 + b A am), m = — + mh (2.8) 

The axial current renormalization constant Za has been computed non perturbatively for the 
0{a) improved theory in [6]. For what concerns the improvement coefficient b A , we know from 
ref. [7] that already at (3 = 6.4 the discrepancy between the one loop calculation and a non- 
perturbative one is of the order of few percent. For this reason we assume for 6^ the perturbative 
value quoted in [8] introducing, in our final result, a systematic error that is below 1%. For large 
times, eq. (2.5) could be used in order to determine Mhe, while the decay constant fhe could be 
extracted from the ratio 



Z A Up)_ ^ l /M( M ML 3 ) l/2 e -(,.-T/ 2 )M u 

Ti 2 



1 + riA*e- XoA + rf A e- {T - Xo)mG 



(2.9) 



In a finite size time extension, the asymptotic expansions (2.5-2.6) are in general not valid, and 
it's impossible to disentangle the lowest state contribution to the correlation functions from the 
excitations due to the higher states. We define masses and decay constants at a value of xq 
which is a fixed fraction of the total time extent of the lattice, e.g. xq = T/2. This affects the 



finite size observables with spurious contributions from the excited states, if T is not sufficiently 
large, but the recursive procedure which connects the small volumes to the large ones through 
the step scaling function, also connects small times to large times, and the final result comes 
out to be projected onto the fundamental state. For this reason we choose to measure masses 
and decay constants at x = L = T/2, through the equation 



All the simulations are done in the quenched theory, and the connection to physical units is 
done fixing the scale ro [9] at r = 0.5 fm. This is a convenient choice for quenched QCD, where 
the ratio r /a has been computed for a wide range of /3's with high precision [10,11] and can 
be connected at higher values of (3 with the behaviour expected from asymptotic freedom (see 
ref. [19] for details). 

2.2 Heavy-light mesons on a finite volume 

An important aspect of the calculation is the tag of the heavy-light meson states at finite 
volume. This is required for both the computation of the decay constant fhe(Lo) on the smallest 
volume and for the evolution steps to the larger ones. The identification has to be done in 
terms of a physical quantity which is independent from the volume and we choose it to be the 
Renormalization Group Invariant (RGI) quark mass. The recipe we propose is the following: 
first of all, we monitor the quark masses via the axial Ward identity 



Then, we connect the Ward Identity masses mwi to the renormalization group invariant (RGI) 
masses M through the relation 



fht{L) = 



2£a U{L) 



(2.10) 



d^AUx) = (m e WI + m^j)P(x), 



(2.11) 



M = Z M (g ) m WI (g ) = Z M (g ) 



(1 + b A am q ) 
(1 + bparriq) 



mwi{go) 



(2.12) 



where 



ZM{go) 



M Z A {g ) 



1 



(2.13) 



m(n) Z P (g ,n) 
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The renormalization constants Zp and Zm have been determined in [12]. Regarding to the 
improvement coefficients 6^ and bp, for a sufficiently small subtracted quark mass in lattice 
units am q , what is really needed is the difference 6 a — bp, which is known non perturbatively 
from [13]. The identification of the heavy-light meson states on a finite volume is done by 
expressing all the observables as functions of the RGI masses and extracting their values at the 
physical points computed in literature and reported in Tab. 1. 

Table 1. RGI masses 

/ M f (GeV) ref. 

s 0.138(6) [14] 

c 1.684(64) [15] 

b 7.01(3)(10) [16] 



2.3 Convergence of the inversion algorithm for the heavy quark propagator 

Propagators are computed using the BiCGStab algorithm with SSOR preconditioner for the 
inversion of the Dirac- Wilson operator plus the clover term [17]. The inversion is done with 
32bit arithmetic both for the light and the heavy quark, and in principle one could ask weather 
rounding effects are present for the heavy quark case, where the exponential decay of the propa- 
gator is quite steep. As it will be explained in sec. 3, the parameters of all the simulations have 
been chosen so to keep an upper bound on the bare heavy quark mass such that omj, < 1/3. In 
this situation the inversion is expected to be safe. Nevertheless, the quality of the inversion can 
be monitored with two simple checks. The first one follows from the observation that the ex- 
ponential decay of a heavy quark propagator in a non trivial fixed gauge background has small 
fluctuations around the tree-level path due to the heaviness of the quark, and the rounding 
effects can be reliably monitored with an analytic comparison at tree-level. The second one 
is suggested by the fact that the quarkonia states hh can be accomodeted on a small volume 
without sensible finite size effects, and the rounding on the propagator can be monitored looking 
at the mass spectrum for these heavy-heavy states. All the simulations we made did pass the 
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two checks. 



3. The specific calculation and the results 

The results are obtained at finite lattice spacing. The size of the smallest volume follows from 
the decision of making our estimate for the finite size Jb on a 48 x 24 3 lattice with a cutoff of 
about clq 1 ~ 12 GeV. The value of the bare coupling for this lattice spacing has been obtained 
from a fit in ref. [19]. The procedure fixes /3(ao) = 7.3 and the physical volume L = 0.4 fm. 
On this lattice, we simulate heavy quark masses up to 0.3 in lattice units, corresponding to bare 
physical masses slightly above 4 GeV. Indeed, as a general caution against large lattice artifacts, 
at all (3 values we take the maximum heavy quark mass in lattice units of the order of 0.3. 
The first E (we distinguish between the continuum step function a and the one at finite lattice 
spacing E) goes from the volume of 0.4 fm to the one of 0.8 fm. In terms of lattice points, we 
go from 12 to 24, and we have to match the starting volume of 0.4 fm with a resolution which is 
half of the one used for a correct estimate of the bottom quark propagation. According to our 
caveat, it follows that the maximum bare quark mass that we can achieve is correspondingly 
halved, i.e. of about a couple of GeV at a bare coupling (3 = 6.737. We make a further iteration 
with a second E going from 0.8 fm to 1.6 fm, where our investigation of heavy quark masses 
stops at the order of the charm quark mass. The corresponding bare coupling is (3 = 6.211. The 
finite volume effects for this second evolution step are small enough to make the neglection of 
the residual volume effects a safe assumption, that however can be tested explicitly. 

The plots in Figs. 1 and 2 show the dependence of E upon the heavy RGI quark mass M RGI for 
the two volume jumps and provide evidence for a plateau of insensitivity to heavy quark masses: 
the first three sets of data represent the measured values of E at fixed values of the light quark 
mass M RGI , and the other two have been obtained from a linear extrapolation in M RGI to the 
down and strange RGI quark masses reported in Tab. 1. The detail plots show a fit, in the 
region of large quark masses only, to the M RGI dependence of E reported here against 1/M RGI 
and confirm the validity of the expansion, given the small slope of the 1/M RGI correction. The 
dependence upon the light quark mass for fixed heavy quark masses is given in Figs. 3 and 4. 
These figures are the core of this paper and support the procedure proposed. 
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Fig. 1. Step scaling function So. 4-0. 8 for the evolution of fhi from 0.4 fm to 0.8 fm at (3 = 6.737. 
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The statistical errors are computed by a jacknife method. More specifically, the errors for the 
values of E at the physical points, i.e. at the values of the charm and of the bottom RGI quark 
masses, are obtained by making two independent jacknife fits to the 1/M^ GI dependence to the 
numerator and to the denominator of eq. (1.1). In this way, we avoid dealing with results at 
various mass values correlated by the same set of background gauge configurations. The final 
error on the E's is obtained by combining quadratically the relative statistical errors resulting 
from the two jacknife fits. 

The finite size value of fs is obtained by a calculation on the highest resolution lattice. The 
RGI bottom quark mass, according to the previous section, is obtained from the equation 

M = Z M (g ) m WI (g ) (3.1) 

In order to obtain the renormalisation constant Zm(9o) at /3 = 7.3 and (3 = 6.737, we have used 
a safe interpolation of the pseudoscalar renormalisation constant Zp (g , fi) at a value of /j, three 
times the reference value used in eq. (6.8) of ref. [12]. The value for fs that we obtain is 

f B ( 0.4 fm ) = 483(4) MeV (3.2) 
By using the values of E for the b quark at constant RGI mass marked in figures 3 and 4 

££4-0.8 = 0-401(4), = 0.88(4) 

we obtain our estimate of fs on the large volume: 

f P B hys = Sb{ 0.4 fm ) • E&_ . 8 • S&^.e = 170(11) MeV (3.4) 
where the error quoted in the previous equation is statistical only. 

The systematic errors can be partly ascribed to specific approximations used in the present 
computation that can be eventually removed, and partly to the uncertainty in the extrapolation 
in the heavy quark mass of finite size effects, inherent to the method proposed. 
To the first class belong the errors related to a finite lattice spacing both for the step scaling 
function and for the finite size decay constant. The former introduce a new dimensionful variable 
aL in the finite lattice spacing step scaling function. The use of an non-perturbatively improved 
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(3.3) 



action, but for surface counterterms that are evaluated in perturbation theory, makes the error 
at most of 0(a lattlce a/L) or 0(a lattlce ciAqcd) i-e. of a few percent. Notice that lattice artefacts 
related to the heavy quark mass alone cancel in the ratio defining the step scaling function at 
diffenet lattice sizes, but at the same values of the cutoff and of the quark mass. 
Preliminary results from a calculation where the continuum limit is estimated indicate that such 
is the case [18]. 

We estimate an uncertainty for each finite lattice spacing step scaling function of about 2%. 
The lattice artifacts of order (aMh) 2 remain in the determination of the finite size With 
our restriction on the maximum value of the heavy quark mass in lattice units, we can limit 
this uncertainty to less than 10%. The overall effect of finite lattice spacing on the quantity in 
eq. (1.2) is not expected to exceed 12 — 13%. 

A second source of uncertainty derives from our estimate of the lattice spacing at large (3 obtained 
from the asymptotic freedom fit of ref. [19]. This produces a variation of fs on the small volume 
in an indirect way. If the lattice spacing, say, is larger than estimated, the value of fs translated 
in physical units is accordingly smaller. However, in such a case, the volume used is larger than 
the expected 0.4 fm: one must repeat the calculation at higher beta, on the matched physical 
volume, where the value of /# is higher, because of finite size effects. The variation induced 
by an error in the lattice spacing depends upon the finite size dynamics. In order to estimate 
it, one has to get the error on (3 for a fixed lattice spacing and perform test simulations within 
the error range. From ref. [19] this amounts to a 0.2% error. We have performed simulations at 
(3 = 7.2 and of (3 = 7.3, a variation range ten times bigger than the error quoted, and obtained a 
variation of /b on the small volume much below our statistical error that must then considered 
a generous upper bound on this effect. An independent estimate could come from bottomonium 
spectroscopy on a finite volume that would also supply an addidtional estimate of the bottom 
quark mass. 

A minor source of uncertainty, negligible and anyway removable, come from our derivation of 
the renormalisation constant that determines the RGI invariant mass from the ward identity 
mass at the highest (3. 

To the first class finally belong the residual finite volume effects beyond the lattice size of 1.6 fm 
that we have simulated. This is a volume considered safe for numerical simulations of light 
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quark spectroscopy and we do not expect residual corrections. As already mentioned, a specific 
test can anyway be made by calculating the next step function. 

To the second class of errors belong the ones deriving from the extrapolation of the step function 
to values of the heavy quark mass higher than the ones simulated, i.e. the validity of the 
asymptotic expansion of eq. (1.2). This can be partly eliminated by running more quark mass 
values and by constraining further the fit. Our data lign on a straight line very well. We estimate 
an error from a parabolic fit through our three points of about 1%. 

The overall error on the number fs coming from the removable systematic uncertainties is of 
about 13% and of at most 2 — 3% from the ones deriving from the unavoidable extrapolation in 
the heavy quark mass, leading to a global uncertainty of about 25 MeV of which about 20 are 
removable while 5 stay with the method: 

f phys = 170 ( 11 )( 5 )(22) MeV (3.5) 

As a check of the whole procedure, we have made the estimate for the charm quark case, whose 
RGI mass from ref. [15] provides a very good fit to the spectroscopy on the largest volume at 
13 = 6.211: 

M D = 1.814(6) GeV, M Vc = 2.881(2) GeV (3.6) 

For the charm quark, the value of the decay constant fu coming from the finite size procedure 
can be compared with the value obtained directly on the large volume L = 1.6 fm at (3 = 6.211. 
The finite size decay constant is 

f D ( 0.4 fm ) = 634(6) MeV (3.7) 

while the step scaling functions are 

Sol-o.8 = 0.397(4), Z c d 8 _ ls = 0.81(2) (3.8) 

The comparison between the two determinations is 

f D hVS = fn( 0.4 fm ) • Eg d 4 _ . 8 • Eg* _ 1-6 = 204(9) MeV 

(3.9) 

f D ( 1.6 fm ) /3=6 .2ii = 208(6) MeV 
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The agreement of the two numbers gives us confidence on the reliability of the fs result. We 
have also extracted the value of fs s by the same procedure, with the finite size decay constant 
and the step scaling functions for the two jumps given by 

f Bs ( 0.4 fm) = 490(4) MeV, s£ s 4 _ . 8 = 0.417(3), E b %_ ls = 0.94(3) (3.10) 

The infinite volume result for this decay constant is 

f phys = 192 (9)(5)(24) MeV (3.11) 

The light quark mass dependence of the volume effects is larger when the volume is large enough 
to resolve the difference between a strange and light quarks. 
Other results coming from our calculations are: 

f V B a ys /fB ys = 1.13(2)(1) (3.12) 

and 

fgv'/fgw = 1.10(1)(1) (3.13) 

where the first error is statistical and the second comes from the uncertainity in the light quark 
extrapolations. 

Finally, we have explored the validity of the phenomenological ansatz of eq. (1.3) for E and 
Fig. 7 shows Eq.4-0.8 as a function of the reduced mass. A reasonable scaling is obtained setting 
Hi = rrii + M , where M = 0.5 GeV, not far from Aqcd- 

We stress again that an estimate of the bottom quark mass and of the lattice spacing at the 
smallest volume can also be made independently by a fit to the bottomonium spectroscopy on 
the smallest size, highest resolution lattice, under the hypothesis that such flavourless states 
have a much smaller radius than the flavoured mesons and do not suffer from finite size effects. 
The difficulty of disentangling the contribution of excited states that, measured on the scale of 
the heavy state, are almost degenerate in mass, could be overcome by optimizing the SF source 
to obtain the best projection on the fundamental state. 

The major systematic effects present in this exploratory calculation can be eliminated, while the 
unavoidable errors due to the extrapolations in the heavy quark mass are strongly suppressed 
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by the manifest insensitivity of finite size effects to the heavy quark mass, which is the main 
result of this paper. All the steps of the calculation always deal with physical quantities properly 
renormalized in massless lattice QCD. 

The method proposed can be generalized to problems characterised by two very different mass 
scales, if the decoupling of the large mass scale from the low scales of non-perturbative QCD 
dynamics holds true. This appears to be the case in the example discussed and is somehow 
supported by the wide success of the predictions of perturbative QCD calculations for hard 
processes that are insensitive to the dressing mechanism of quarks and gluons into standard 
hadronic final states. 
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